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The  various  steps  of  the  research  are  as  follows  : 

-  Development  of  a  picosecond  source  providing  tunable  pulses 
(570-700nm)  of  1,5  mJ  energy. 

-  Construction  of  cells  (4  photon  resonance  in  toercury)  and 
heat  pipes  (2  photon  resonance  in  lithium,  single  and  double  heat  pipes,  and 
cyllndrlc  heat  pipes  for  ionisation  measurements). 

-  Development  of  a  data  acquisition  system  capable  of  averaging 
out  data  as  a  function  of  relative  phase,  delay  and  amplitude  of  the  exciting 
pulse  sequence.  The  averaging  operations  have  to  be  performed  at  a  data  rate 
of  up  to  20  Hz. 

-  Measurement  of  transmission  (tmiltiphofon  absorption)  for  a  two 
pulse  sequence,  as  a  function  of  the  delay  and  phase  difference  between  the 
two  pulses,  and  of  the  pulse  energy.  The  measurements  are  to  be  conducted 

in  dye  vapors,  lithium  vapor  and  mercury  vapor  (cases  corresponding  respec¬ 
tively  to  a  single, double  and  four  photon  resonance  ). 

-  Measurements  of  third  and  fifth  harmonic  generation  in 
conditions  of  optimal  multiphoton  excitation  determined  in  the  previous 


measurements . 


progress , 


resonances 


-  Simultaneous  theoretical  analysis  of  the  measurements  in 


-  Implementation  of  the  concept  of  enhancement  by  "induced 


III.  STATUS  OF  THE  RESEARCH 


III. I  The  source  and  the  Data  Requisition  system 


The  development  of  the  source  has  been  completed,  providing 
tunable  picosecond  pulses  of  1,5  mJ  energy,  at  a  repetition  rate  of  up  to 
20  Hz.  The  laser  and  associated  data  acquisition  system  developed  for  the 
experiments  of  coherent  multiphoton  propagation  are  described  in  a  publi¬ 
cation  (submitted  to  Applied  Optics,  Included  as  Appendix  A).  First  we 
summarize  the  performances  of  the  system,  as  described  in  the  latter 
publication.  Next,  we  will  focus  on  some  shortcomings,  and  describe  the 
solution  that  are  being  implemented. 


III. 1.1  Amplifier  characteristics 


We  use  a  three  stage  amplifier  (two  pumped  tranversely,  one 
longitudinally).  Coupling  between  the  first  stages  is  minimized  by  synchro¬ 
nizing  the  first  stage  at  the  leading  edge  while  the  second  stage  is  synchro¬ 
nized  at  the  trailing  edge  of  the  pump  pulse.  The  amplified  spontaneous 
emission  (ASE)  after  the  last  stage  is  approximately  1%  of  the  pulse  energy, 
without  the  use  of  any  other  "isolator"  than  a  direct  vision  prism  between 
the  two  first  stages. 

plified  beam,  which  is  essential  to  perform  the  critical  alignement  of  the 
interferometers  used  for  the  analysis  of  time  and  space  coherences. 

Beam  Profile  Careful  positioning  of  the  pump  beams  relative  to  the 

amplified  beam  in  the  last  two  stages  provides  complete 
control  over  the  spatial  intensity  distribution,  from 
"square"  to  "Gaussian". 

Spatial  Coherence  of  the  amplified  pulse  is  demonstrated  by  achieving 

uniform  extinction  of  the  beam  in  a  Mack  Zehnder  inter¬ 
ferometer. 

Temporal  Coherence  The  measurement  of  temporal  coherence  is  crucial  to  the 

experiments  of  coherent  propagation,  and  has  been  used  to 
test  the  laser  sources,  the  scanning  technique,  and  the 
data  acquisition  system.  The  measurements  detained  in 
Appendix  A  (ref.l  and  2)  indicate  that  the  amplifier  does 
not  introduce  any  phase  modulation.  The  amplified  pulse 
can  best  be  represented  by  a  Gaussian  temporal  profile 
with  a  phase  modulation  <(>(t)  proportional  to  the  pulse 
intensity  (<)>»1.43  E^).  Such  a  determination  cannot  be  made 
with  conventional  spectral  measurements,  since  such  a 
chirp  corresponds  to  an  increase  in  the  bandwidth  duration 
product  of  only  10%. 

III. 1.2  Data  Acquisition 

A  Mack  Zehnder  interferometer  is  used  to  create  the  sequence  of 
two  pulses  with  continuously  increasing  relative  phase  and  delay  (as  needed 
for  the  experiments  of  coherent  propagation,  as  well  as  for  the  analysis 


of  pulse  coherence).  The  data  acquisition  system  is  needed  to  resolve  two 
conflicting  requirements  : 


a)  the  scanning  speed  has  to  be  fast  enough  to  ensure  continuous 
frictionless  motion  of  the  delay  arm  (minimum  speed  :  one  wavelength  / 2s)  ; 

b)  averaging  of  a  large  number  of  data  points  per  wavelength 
increment  of  the  delay  arm,  which  would  require  even  slower  scanning  speeds. 

The  solution  that  has  been  first  implemented  (described  in 
Appendix  A  pg  A8-A10  ref. 1,2)  is  to  increase  the  delay  at  a  regular  rate  of 
0.l6pm/s,  to  divide  each  fringe  by  an  equal  number  of  phase  intervals,  and 
to  average  data  of  corresponding  phase  over  a  set  of  —for  instance- 
100  fringes. 

Despite  its  shortcomings  discussed  below,  this  method  enabled  us  to  obtain 
accurate  and  simultaneous  measurements  of  Intensity  and  interferometric 
autocorrelations,  and  to  perform  as  well  nonlinear  transmission  measurements 
through  lithium  vapor. 

III. 1.3  Latest  improvments 

It  was  established  that,  while  delay  increments  of  1000  pro 
occured  at  equal  time  intervals  (of  lh,44min, 10s) ,  the  translation  speed 
showed  unacceptable  fluctuations  over  distances  of  a  few  microns.  These 
speed  fluctuations  cause  the  data  in  successive  periods  to  be  assigned  to 
the  wrong  "phase  channel",  resulting  in  a  decrease  of  the  "phase  contrast" 
(between  data  corresponding  to  different  "phase  settings")  to  be  observed. 

Inchworm  An  "inchworm"  was  purchased  from  "Burleigh  Instruments",  with 

matching  funds  from  the  University.  According  to  the  manufactor's 
advertisements, this  piezoelectrically  driven  device  should  enable 
us  to  reproduclbly  scan  a  delay  with  increments  of  only  0.01  pm. 
Our  measurements  instead  showed  a  nonlinearity  such  that  the 
number  of  steps/period  changed  more  than  50%  over  only  5  periods. 
This  approach  had  to  be  abandonned. 

Simultaneous  measurement  of  the  distance 

The  problems  that  arose  from  the  previous  approach  originate  from 
the  use  of  time  rather  than  space  as  the  primary  variable. 


The  interference  maxima  (or  "in  phase  sequences")  did  not  occur 

at  equal  time  intervals  because  of  speed  fluctuations.  This 

problem  is  eliminated  if  the  absolute  distance  is  continuously 

measured  with  subwavelength  accuracy.  This  is  done  by  sending 

a  He-Ne  laser  beam  through  the  same  optical  delay  (Mack  Zehnder 

Interferometer)  as  the  picosecond  pulses,  and  monitoring  the 

intensity  of  the  transmitted  beam  to  determine  the  distance. 

In  order  to  minimize  the  computing  time,  conversion  of  HeNe 

fringe  intensity  to  distance  was  first  made  by  comparing  the 
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measurements  to  a  table  of  cos(-y  Z^)  values.  The  purchase  of 
a  new  fast  arithmetic  card  for  the  computer  enabled  us  to  define 
the  distance  by  using  directly  an  arc  cos(I)  subroutine.  The 
computed  distance  is  used  to  determine,  in  each  successive  wave¬ 
length  increment,  the  exact  "phase  channel"  to  which  a  particular 
data  point  belongs.  The  data  are  next  averaged  per  channel,  for 
a  large  number  of  wavelengths,  as  described  in  Appendix  A  (ref.l 
and  2)  . 

III. 2  Heat  Pipe  and  Cell  Construction 

The  rather  slow  scanning  process  (circa  15  hours)  imposes  severe 
contralnts  not  only  on  the  amplifier  chain  stability  but  also  on  the  heat  pipe. 
Obviously  the  heat  pipe  has  to  be  perfectly  leak  tight  but  in  addition  the 
outgasslng  of  the  walls  has  to  be  continuously  and  selectively  (l.e.  without 
pumping  the  He  buffer  gaz)  eliminated  during  that  15  h  period.  Another  problem 
is  a  condensation  of  solid  lithium  at  the  edges  of  the  hot  zone.  The  resulting 
accumulation  over  several  hours  of  an  annular  ridge  of  solid  was  sufficient 
to  obturate  the  optical  path  of  the  original  heat  pipe.  To  solve  this  problem 
we  constructed  another  (single)  heat  pipe  of  2"  inside  diameter  (instead  of 
1"  for  the  original  one).  Two  photon  propagation  measurements  (appendix  B) 
could  be  performed  with  this  heat  pipe  over  extended  periods  of  time.  The 
construction  of  another  double  heat  pipe  (of  2"  inner  diameter)  is  about  to 
be  completed,  to  be  used  for  the  measurements  of  harmonic  generation  in 
mixtures  of  lithium  and  magnesium. 


In  order  to  perform  comparative  measurements  of  coherent  propa¬ 
gation  and  harmonic  generation  in  conditions  of  single  and  four  photon  reso¬ 
nances  .vacuum  cells  have  been  constructed  for  dye  vapors  (malachite  green) 


and  mercury  vapor  (for  4  photon  transition).  Because  of  the  relatively  high 
vapor  pressure  of  these  samples,  simple  pyrex  cells  have  been  used,  with 
fused  silica  or  LiF  windows.  The  sealed  cells  are  uniformly  heated,  with  a 
slightly  higher  heating  power  applied  to  the  windows,  to  prevent  condensation. 

In  order  to  study  the  ionisation  losses  associated  with  two  photon 
excitation  of  lithium,  a  cylindrical  heat  pipe  has  been  designed  and  is 
presently  being  assembled.  The  upper  and  lower  halves  of  this  heat  pipe  are 
electrically  Isolated  making  them  useable  as  electrodes. 

HI. 3  Theory 

The  computer  controlled  data  acquisition  data  acquisition  system 
provides  a  continuous  recording  of  absorption-versus  phase  and  delay.  In  order 

to  match  this  large  set  of  data  with  the  theoretical  model,  a  computer  program 
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has  been  made  that  simulates  the  complete  experiment  (Appendix  B) .  In  order 
to  bring  the  computing  time  under  reasonable  limits,  the  propagation  problem 
has  been  linearized  (thin  sample  approximation) . 

A  theoretical  study  of  multiphoton  excitation  under  sequential 

excitation  by  coherent  pulses  is  detailled  in  a  publication  submitted  to  the 
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Journal  of  Chemistry  and  Physics  (Appendix  C) .  It  Is  shown  that,  by  using 
a  sequence  of  Identical  pulses  with  well  defined  phase  relationship,  one 
combines  the  spectral  selectivity  offered  by  a  regular  pulse  train,  with  the 
high  peak  intensity  of  the  short  pulses,  which  are  needed  to  excite  higher 
order  processes.  The  case  of  random  pulse  sequences  is  also  considered,  and 
shown  to  offer  an  accurate  methqd  to  study  the  influence  of  source  coherence 
on  a  multiphoton  process. 
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ABSTRACT 

We  describe  a  high  power  picosecond  dye  laser  system,  providing  tunable 
1.5  mJ  pulses,  with  a  uniform  wavefront  and  Gaussian  spatial  profile.  With  the 
aid  of  a  computer  controlled  data  acquisition  system,  we  have  analyzed  for  the 
first  time  the  temporal  coherence  properties  of  such  a  system. 


I.  INTRODUCTION 

We  report  the  development  of  a  picosecond  pulse  dye  laser  osci.lator- 
amplifier  system,  designed  for  experiments  of  resonant  coherent  propagation  [i]. 
Such  experiments  put  stringent  requirements  on  the  properties  of  the  laser 
system,  most  of  which  will  be  reviewed  briefly  in  section  II.  We  have  been  able 
to  meet  these  requi rements,  using  a  commercially  available  synchronously  pumped 
dye  laser  in  combi  nati  on  wi  th  a  home  built  ampl  i  f  i  er.  As  we  wi  1  1  expl  ai  n  i  n 
section  III,  this  amplifier  consists  of  three  stages,  being  pumped  by  a  fre¬ 
quency  doubled  Nd:YAG  laser.  The  first  two  stages  are  pumped  transversely, 
while  the  last  stage  is  pumped  longitudinally.  We  succeeded  in  achieving  a  high 
gain,  while  avoiding  a  build-up  of  amplified  spontaneous  emission  (ASE),  without 
the  use  of  spatial  filters  or  saturable  absorbers.  Accurate  synchronization  of 
the  various  pulses,  combined  with  a  careful  optimization  of  the  pump  and  ampli¬ 
fied  beam  geometry,  ensures  that  the  energy  of  the  pump  is  channeled  into  the 
picosecond  pulse  rather  than  into  ASE.  The  advantage  of  our  transverse- 1 ongitu- 
dinal  pumping  geometry,  over  purely  longitudinally  pumped  systems,  is  the  abil¬ 
ity  to  control  the  spatial  beam  profile,  without  compromising  the  output  energy. 
In  section  IY,  we  wi  1 1  present  interferometric  autocorrel ation  measurements 
r2,3],  performed  for  the  first  time  on  such  a  picosecond  pulse  dye  laser  oscil- 
lator-ampl  ifier  system.  These  measurements,  for  which  a  new  kind  of  data  acqui¬ 
sition  system  needed  to  be  developed,  show  that  our  amplifier  does  not  change 
significantly  the  temporal  properties  of  the  dye  laser  pulses.  Preliminary 
measurements  of  coherent  two  photon  absorption  in  lithium  demonstrate  that  this 
system  meets  the  stringent  coherence  requirements  of  a  multiphoton  coherent 
propagation  experiment. 
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II.  GENERAL  DESIGN  CONSIDERATIONS  FOR  A  PICOSECOND  PULSE  DYE  AMPLIFIER 


Even  though  we  designed  our  amplifier  for  achieving  a  particular  goal 
(saturation  of  the  2s-4s  transition  of  lithium  at  571.2  nm),  the  main  specifica¬ 
tions  are  desirable  requirements  for  a  large  number  of  experiments:  (1)  high 
output  energy;  (2)  minimum  ASE;  (3)  good  beam  quality  for  interferometric 
experiments,  i .e.,  uniform  beam  prof i  1  e  and  flat,  di storti on  free  wavefronts; 
(4)  no  temporal  phase  modulation  effects;  (5)  both  short  term  (pulse  to  pulse) 
and  long  term  stability;  (6)  availability  of  a  c.w.  pilot  beam  propagating 
through  the  system,  for  the  ease  of  alignment  of  the  target  experiment.  Each  of 
these  requirements  dictated  one  or  more  basic  design  features  of  the  amplifier 
to  be  bui  It,  as  we  wi  1  1  expl  ain  now.  (See  e.g.,  ref.  4  and  5,  and  additional 
references  therein  for  a  theoretical  analysis.) 

Obviously,  a  high  output  energy  calls  for  an  optimum  use  of  the  pump 
radiation,  which  in  turn  requires  the  minimum  beam  diameter,  in  each  amplifier 
stage,  compatible  with  a  maximum  unsaturated  gain.  Maximizing  the  gain  for  a 
given  input  also  requires  minimal ization  of  ASE.  This  can  be  achieved  through: 
(1)  accurate  timing  of  the  pump  pul se  rel  ati ve  to  the  pi co second  pul se  to  be 
amplified  in  each  stage,  thereby  minimizing  the  overlap  between  the  gain  pulse 
in  one  cell  and  the  ASE  from  the  adjacent  cell;  (2)  accurate  geometrical 
matching  of  the  pump  beam  and  the  picosecond  dye  laser  beam,  in  order  to  maxi¬ 
mize  uniform  pump  depletion  by  the  latter  in  each  cell;  (3)  appropriate  choice 
of  dyes.  The  choice  of  dyes  is  further  limited  by  our  requirement  that  the 
system  should  be  transparent  for  the  unamplified  dye  laser  pulses,  which  serve 
the  purpose  of  a  c.w.  pilot  beam.  This  last  condition  is  particularly  important 
in  applications  which  require  a  critical  alignment  of  the  amplified  beam.  For 
our  applications,  the  dye  laser  beam  has  to  be  sent  first  along  the  axis  of  an 
interferometric  delay  line,  and  subsequently  through  a  long  but  narrow  heat 
pipe,  an  operation  that  requires  a  c.w.  beam.  The  requirement  of  an  aberration 
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free  temporal  and  spatial  beam  profile  imposes  the  following  design  features: 
(1)  short  dye  cells  in  order  to  minimize  pulse  broadening  and  phase  modulation 
effects  due  to  the  solvent  [6,7];  (2)  the  absence  of  spatial  filters  and  satur¬ 
able  absorbers  in  order  to  avoid  wavefront  distortion  in  subsequent  amplifiers; 
(3)  a  choice  of  solvents  to  obtain  maximum  thermal  conductivity;  (4)  careful 
lateral  adjustment  of  the  pump  beam  in  the  longitudinally  pumped  cell,  to  pro¬ 
perly  shape  the  spatial  beam  profile  of  the  amplified  dye  laser  pulses. 
Finally,  stability  requirements  put  stringent  conditions  on  the  pump  source  for 
the  amplifier:  (1)  short  term  stability  requires  a  jitter  free  synchronization 
between  dye  laser  and  amplifier  pump;  (2)  long  term  stability  merely  imposes  a 
constant  thermal  lens  effect  for  the  pump  laser. 


III.  SYSTEM  DESCRIPTIONS 

The  oscillator  is  a  Spectra-Physics  dye  laser,  synchronously  pumped  by  a 
Spectra-Physics  Ar  ion  laser,  which  is  actively  model ocked  at  82  MHz.  For  the 
gain  medium  we  use  Rh  6G  in  Ethylene  Glycol  (concentration  2  x  10"^  M). 
Although  Rh  560  would  yield  a  slightly  higher  gain  at  571  nm,  we  prefer  instead 
to  use  Rh  6G  because  of  its  better  chemical  stability.  (The  Rh  560  solution  has 
to  be  replaced  every  week.)  Wavelength  tuning  around  571  nm  is  achieved  with  a 
3-plate  birefringent  filter.  With  a  30%  output  coupler,  this  system  generates  a 
continuous  train  of  3.6  nj  pulses  with  a  duration  of  6.2  ps  FWHM. 

Our  three  stage  dye  amplifier,  sketched  in  Fig.  1,  is  pumped  by  a  frequency 
doubled,  Q-switched  Nd:YAG  oscillator-amplifier  (Quanta  Ray,  DCR-1).  This  sys¬ 
tem  delivers  280  mJ  pulses  (530  nm)  of  8  ns  FWHM,  at  a  repetition  rate  of  10  Hz. 
The  polarization  of  this  beam  is  orthogonal  to  the  polarization  of  the  dye  laser 
beam.  A  1.5  m  optical  delay  separates  the  dye  laser  oscillator  from  the  first 
amplifying  cell,  in  order  to  prevent  oscillation  of  the  latter  with  the  output 
mirror  of  the  former. 


a .  Beam  Geometry 

To  prevent  ASE  from  building  up  between  stages,  a  direct  vision  prism  is 
inserted  between  the  first  and  second  stage.  This  is  the  sole  "optical  isola¬ 
tor"  used  in  our  system,  in  contrast  to  other  set  ups  which  need  in  addition 
saturable  absorbers  and  pinholes  f4,8,9].  Particularly  important  in  this 
respect  is  the  expanding  beam  geometry  throughout  the  system,  not  only  to 
increase  the  effectiveness  of  the  prism  in  discriminating  against  ASE,  but  also 
to  maximize  the  overall  gain  of  the  amplifier.  Starting  with  a  diameter  of  less 
than  1  mm  in  the  first  cell,  the  beam  diameter  is  increased  3  times  in  the 
second  cell,  and  further  expanded  drastically  in  the  last  stage.  The  latter  is 
inserted  between  two  lenses  of  a  collimating  beam  expander  with  an  aperture 
ratio  of  5.  This  particular  design  configuration  of  the  last  stage  serves  a 
dual  purpose.  The  expansion  of  the  dye  laser  beam  is  required  in  order  to  keep 
the  energy  density  at  571  nm  below  the  saturation  energy  density.  On  the  other 
hand,  the  converging  geometry  of  the  pump  beam  is  needed  to  maintain  a  suf¬ 
ficient  level  of  excitation,  despite  absorption  in  the  dye.  One  should  bear  in 
mi  nd  that  the  performance  1  imi  ts  of  the  dye  ampl  i  f  i  er  are  def  i  ned  i  n  the  last 
stage.  Indeed,  the  saturation  energy  density  over  the  output  cross  section  of 
the  last  stage  is  the  maximum  energy  per  pulse  that  the  system  can  deliver 
(e.g.,  1  mJ/cm^  for  Rh  6G)  [10].  Because  of  its  large  cross  section,  this  last 
stage  requires  most  of  the  energy  (77%)  of  the  pump  laser,  while  yielding  only  a 
low  gain  (factor  8).  Of  the  remaining  pump  energy,  7%  is  used  for  the  first 
stage  (gain:  103),  and  16%  for  the  second  stage  (gain:  50). 

b.  Choice  of  Oyes 

For  the  transversely  pumped  stages,  we  use  1  cm  long  flow  through  cells,  at 
Brewster  incidence.  The  best  dye  candidate  for  a  pump  wavelength  near  530  nm, 
is  Rh  6G  (highest  gain  at  571  nm,  longest  lifetime).  However,  the  absorption 


band  of  Rh  6G  extends  wel  1  beyond  571  nm.  Therefore,  Rh  560  has  to  be  used  in 
the  first  two  stages  (where  the  concentration  is  necessarily  high),  in  order  to 
get  a  reasonably  large  transmission  of  the  unamplified  dye  laser  pulses  through 
the  system.  It  was  determined  empirically  that  the  following  dye  combinations 
in  the  first  and  second  cell  yield  the  best  results  (maximum  pump  depl  etion, 
minimum  ASE,  and  acceptable  loss  at  571  nm):  for  the  first  stage,  a  mixture  of 
Rh  560  and  Rh  6G  in  methanol  (respectively  2  x  10"^  M,  and  10“^  M);  for  the 
second  stage  Rh  560  in  methanol  { 10" ^  M).  Finally,  for  the  last  stage,  which  is 
5  cm  long,  we  use  Rh  6G  to  ful  ly  maximize  the  gain.  The  sol  vent  in  thi  s  stage 
is  a  0.5%  aqueous  solution  of  Ammonix  L.O.,  to  minimize  thermal  effects.  For 
the  small  dye  concentration  in  this  final  stage  (1.5  x  10~5  M),  the  absorption 
of  the  unamplified  pulses  is  negligible. 

c.  Synchronization 

In  Section  II,  we  emphasized  the  need  for  an  accurate  and  jitter  free 
synchronization,  in  each  stage,  between  the  pump  pulse  and  the  dye  laser  pulse 
to  be  amplified.  To  reach  this  goal,  we  use  the  model ocker  drive  of  the  Ar  ion 
laser  as  an  absolute  clock.  Through  direct  digital  division,  this  clock  signal 
is  scaled  down  to  approximately  10  Hz.  Next,  this  signal  is  used  to  interrogate 
the  original  clock  signal  in  a  logical  AND  gate,  thus  providing  a  10  Hz  refer¬ 
ence  signal,  which  is  perfectly  synchronized  with  the  master  clock.  A  detailed 
description  of  the  electronic  circuit  has  been  published  elsewhere  [11].  The 
10  Hz  reference  in  turn  is  used  to  directly  trigger  the  Q-switch,  and  after  an 
appropriate  delay,  also  the  flashlamps  of  the  Nd:YAG  oscillator-amplifier.  The 
resulting  jitter  in  the  synchronization  between  the  pump  pulse  and  the  dye  laser 
pulse  to  be  amplified,  is  of  the  order  of  1  ns. 

To  further  minimize  energy  fluctuations  of  the  amplified  pulse,  due  to  this 
small  remaining  jitter,  we  synchronize  the  first  stage  at  the  leading  edge  of 


the  pump  pulse,  while  the  second  stage  is  synchronized  at  the  trailing  edge. 
This  procedure  further  helps  in  reducing  undesirable  coupling  in  ASE  between  the 
initial  stages. 

d.  Performance  of  the  System 

The  output  energy  of  our  amplifier  is  circa  1.5  mj  per  pulse.  Pulse  to 
pulse  stability  is  better  than  20*.  This  is  a  good  result  for  a  chain  of 
amplifiers,  none  of  which  being  saturated.  Using  fresh  dyes,  the  long  term 
stability  is  also  good  (5%  over  a  period  of  10  hours).  For  an  optimum  synchro¬ 
nization,  the  ratio  of  ASE  energy  to  pulse  energy  after  the  last  stage,  is  in 
the  order  of  a  few  percent  only. 

The  spatial  intensity  distribution  of  the  amplified  beam  has  also  been 
examined,  using  an  optical  multichannel  analyzer  (OMA).  Because  the  first  two 
stages  are  pumped  transversely,  one  does  not  expect  a  uniform  gain  profile 
through  the  whole  chain.  Therefore,  one  may  fear  that  the  amplifier  design  may 
not  provide  a  useful  and  smooth  pulse  spatial  profile.  Our  measurements  instead 
have  shown  that,  by  carefully  positioning  the  pump  beams  relative  to  the  ampli¬ 
fied  beam  in  the  last  two  stages,  we  can  achieve  complete  control  of  the  profile 
of  the  output  beam.  An  example  is  shown  in  Fig.  2.  With  the  OMA  in  the  free 
running  mode,  no  shot  to  shot  variations  of  this  profile  were  detected. 

The  quality  of  the  amplifier  output  was  also  examined  with  a  Mack-Zehnder 
interferometer.  Complete  extinction  of  the  beam  over  its  entire  cross  section 
could  be  achieved,  with  the  interferometer  balanced.  This  demonstrates  the 
absence  of  wavefront  distortion  in  the  amplifier  chain. 

IV.  TEMPORAL  PROPERTIES  OF  THE  AMPLIFIED  PULSES 

Equally  important  for  coherent  propagation  experiments,  are  the  temporal 
properties  of  the  amplified  pulses.  We  infer  these  properties  from  second  order 
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intensity  and  fringe  resolved  autocorrelation  measurements,  using  a  Mack-Zehnder 
interferometer  [2,3].  The  main  difficulty  involved  in  these  measurements,  stems 
from  the  fact  that  the  delay  arm  of  the  interferometer  has  to  be  moved  uni¬ 
formly,  with  absolute  subwavelength  accuracy,  over  a  sufficiently  long  distance. 
Indeed,  the  requirement  that  the  motion  of  the  delay  arm  be  continuous,  without 
the  di sconti nui ties  that  would  result  from  the  static  friction  coefficient 
(being  larger  than  the  dynamic  one),  implies  a  minimum  scanning  speed.  On  the 
other  hand,  however,  an  upper  limit  for  this  scanning  speed  is  set  by  the 
repetition  rate  of  the  amplifier  (10  Hz),  and  the  requirement  that  there  be  at 
least  twenty  to  thirty  data  points  per  period,  in  order  to  resolve  individual 
fringes.  Obviously,  both  requirements  are  conflicting.  Moreover,  one  has  to 
realize  that  the  remaining  pulse  to  pulse  instabilities,  will  cause  large  fluc¬ 
tuations  on  top  of  second  order  fringes,  thereby  washing  out  the  fringe  con¬ 
trast.  Averagi ng  out  the  data  over  a  1  arge  number  of  laser  shots,  wi  1 1  only 
further  reduce  this  contrast.  To  solve  these  problems,  we  had  to  develop  an 
adequate  data  acquisition  system. 

a.  Oata  Acquisition  System 

The  hardware  of  our  data  acquisition  system  is  schematically  represented  in 
Fig.  3.  About  1 %  of  the  amplifier  output  is  used  as  a  reference  beam,  the 
remaining  part  is  sent  through  the  interferometer.  Both  the  reference  beam  and 
one  of  the  tput  beams  of  the  interferometer  are  frequency  doubled  and  digi¬ 
tized  if  •"  ^ction  by  photoconducti  ve  p-i-n  diodes  (UDT  100L).  The  other 
output  beam  ,e  interferometer  is  used  for  two  photon  coherent  propagation 
experiments  in  lithium  vapor.  The  electronic  discriminator  in  the  reference 
channel,  rejects  pulses  with  either  too  low  or  too  high  U.V.  energy.  This 
feature  is  built  in,  mainly  in  view  of  our  target  experiment,  which  depends 
critically  on  the  pulse  energy. 


The  forementioned  problem  of  smoothing  the  interference  fringes  is  solved 
in  the  software  of  our  data  acquisition  system.  First,  errors  due  to  amplitude 
fluctuations  of  the  pulses,  are  eliminated  by  computing  the  ratio  of  the  signals 
i n  channel  2  and  1.  Next,  whi 1 e  vary i ng  the  delay  in  the  autocorrel a tor  at  a 
regular  speed  of  0.16  ym/s,  the  resulting  fringes  are  partitioned  over  a  large 
number  of  sets  (each  set  has  typical  ly  100  fringes).  Per  set,  two  different 
averaging  processes  are  performed.  On  one  hand,  all  data  in  each  set  are 
averaged  out,  thereby  washing  out  the  fringes  completely.  For  a  large  number  of 
sets,  this  will  yield  the  conventional  3:1  intensity  autocorrelation  curve.  On 
the  other  hand,  an  average  fringe  per  set  is  constructed,  by  averaging  out  the 
data  in  a  single  periodicity  over  the  whole  set.  An  example  is  shown  in  Fig.  4. 
This  averaging  process  permits  one  to  keep  track  of  small  phase  shifts  over 
consecutive  sets  of  fringes.  For  a  large  number  of  sets,  it  yields  the  fringe 
resolved  8:1  autocorrelation  curve.  In  this  way,  one  thus  obtains  during  a 
single  experiment,  both  autocorrelation  curves. 

b.  Data  Reduction 

By  fitting  the  experimental  data  for  both  the  intensity  and  the  fringe 
resolved  autocorrelation  curves,  it  is  possible  to  determine  accurately  the 
temporal  properties  of  the  amplified  pulses.  From  the  experimental  set  up 
described  above,  we  obtain  two  recordings  of  second  harmonic  signal  versus 
delay:  the  intensity  and  the  interferometric  autocorrelation.  The  intensity 
autocorrel ation  curve  provides  a  first  estimate  of  a  reasonable  pulse  shape  and 
pulse  duration,  which  is  used  as  a  starting  point  in  attempting  to  fit  the 
interferometric  autocorrel ation  measurements.  For  a  given  trial  function  for 
the  complex  electric  field  envelope 


a  numerical  computation  is  made  of  the  interferometric  second  order  autocor¬ 
relation 
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of  which  the  intensity  autocorrelation  I(t)  is  the  subpart, excl uding  the  rapidly 
varying  terms  in  cos  (2*vt)  and  cos  (4*vt)  of  (2): 
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The  computations  are  organized  as  follows:  first,  at  x  =  0,  and  at  t  *  l/2v  the 
first  fringe  maximum  and  minimum  is  computed;  subsequently,  x  is  incremented  by 
at  =  1/v,  to  generate  the  next  pair  of  points,  and  so  on.  In  order  to  improve 
the  accuracy,  an  additional  loop  generates  near  each  fringe  maximum  three  points 
of  F(t  ) ,  not  more  than  two  of  them  being  located  at  the  same  side  of  that 
maximum.  With  a  parabol  ic  fit  through  these  points,  the  computer  determines: 
(1)  the  true  location  of  that  maximum  (Tmax),  (2)  the  correct  at,  which  will  be 
used  to  initiate  the  next  main  loop,  (3)  the  location  of  the  fringe  minimum 
(*min),  and  finally  (4)  the  corresponding  values  of  F(Tma)t),  F(Tm1n)  and 
I ( Tma x ) -  Once  I(t  )  and  both  envelopes  of  F(x)  are  computed  in  this  way,  the 
pulse  parameters  (shape,  asyimaetry,  width,  and  chirp)  arc  adjusted,  until  a  good 
fit  with  the  experimental  data  can  be  reached.  Final  ly,  for  the  best  fitting 
pulse,  the  deviation  of  the  average  frequency 
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and  the  pulse  bandwidth 


are  obtained. 


c.  Results:  Pulse  Characteri sties  Before  and  After  Amplification 

Fig.  5a  shows  the  intensity  autocorrel ation  curve  for  the  amplified  pulses 
(dashed  line).  The  solid  line  represents  a  numerical  fit,  using  a  Gaussian 
pulse  shape  with  6.2  ps  FWHM  (i.e.,  5.27  ps  HW( l/e^)M).  Notice  that  an  inten¬ 
sity  autocorrel  ati  on  cannot  effectively  distinguish  between  a  Gaussian  and  a 
Sech  pul se  shape. 

From  Eq.  (3),  it  is  clear  that  the  intensity  autocorrel ation  is  phase 
independent,  and  therefore  does  not  carry  any  information  about  the  pulse  chirp. 
On  the  other  hand,  as  shown  in  ref.  [2,3],  the  interferometric  autocorrelation 
is  sensitive  to  the  pulse  chirp  f(t).  The  experimental  resul  ts,  obtained  for 
the  1  atter  with  our  ampl  i fied  pul  ses,  are  shown  in  Fig.  5b  (dashed  1  ine).  The 
solid  line  in  Fig.  5b  represents  again  the  best  theoretical  fit:  Gaussian 
pulses  (6.2  ps  FWHM)  with  a  phase  modulation  proportional  to  the  pulse  inten¬ 
sity  ( <>  =  1.43  E^).  For  comparison,  the  8:1  autocorrelation  curve  for  the  same 
Gaussian  pulses,  but  having  instead  a  linear  phase  modulation  (best  fit: 
<l>  =  1.5  (t/5.27)^  is  also  shown  in  Fig.  5b  (dotted  line).  As  a  result,  the 
amplified  pulses  are  not  completely  bandwidth  limited;  according  to  Eq.  (5),  the 
bandwidth  is  79  GHz. 

The  data  shown  in  Fig.  5b  are  much  more  sensitive  to  phase  modulation 
effects  than,  e.g.,  spectral  data.  The  bandwi dth'duration  product  for  the 
nonlinear  chirped  pulses  in  our  case  is  0.49,  or  only  10%  larger  than  the 
minimum  of  0.44  for  unchirped  Gaussian  pul  ses.  On  the  other  hand,  the  FWHM  of 
the  upper  envelope  of  the  interferometric  autocorrelation  curve  for  our  pulses 
is  7.2  ps,  compared  to  9.5  ps  for  unchirped  pulses;  the  lower  envelope  is  even 
more  sensitive  to  phase  modulation  effects:  its  FWHM  changes  by  almost  50%  from 
our  chirped  pulses  (9.5  ps)  to  the  case  of  unchirped  pulses  (14.7  ps). 


The  kind  of  chirp  we  have  detected  for  the  amplified  pulses,  i.e.,  a  chirp 
proportional  to  the  time  derivation  of  the  pulse  intensity,  is  typical  for  self¬ 
phase  modulation  of  the  pulses.  An  interferometric  autocorre 1 ation  of  the 
unamplified  dye  laser  pulses,  as  shown  in  Fig.  6,  clearly  indicates  that  the 
origin  of  the  chirp  is  in  the  laser  oscillator,  rather  than  in  the  ampl  i  f  ier 
chain.  As  a  matter  of  fact,  partial  chirp  compensation  is  taking  place  in  the 
amplifier:  the  interferometric  autocorrel ation  curve  for  the  unamplified  pulses 
reveals  a  phase  modulation  which  is  proportional  to  1.84  E2  (rather  than  1.43  E^ 
for  the  amplified  pulses). 

To  conclude,  our  Interferometric  autocorrelation  measurements,  such  as 
shown  in  Fig.  5b,  have  demonstrated  that  the  phase  distortion  introduced  by  the 
amplifier  is  small  in  the  time  as  well  as  in  the  space  domain.  A  peak  to 
background  of  8:1  can  only  be  obtained  with  spatially  coherent  pulses.  Indeed, 
at  zero  delay,  the  maximum  signal  of  eight  times  the  background,  requires  that 
the  two  halves  of  the  pulse  interfere  constructively  over  their  entire  cross 
section. 


Y.  CONCLUSIONS 

In  summary,  we  have  devel  oped  a  picosecond  pul  se  dye  1  aser  osci  1  1  ator- 
amplifier  system,  tunable  near  570  nm.  At  a  repetition  rate  of  10  Hz,  the 
output  consists  of  1.5  mJ  pulses.  Both  the  short  term  and  long  term  stability 
Is  excel  lent.  Contrary  to  other  simi  1  ar  systems,  we  have  avoided  the  use  of 
either  pinholes  or  saturable  absorbers.  Therefore,  our  system  requires  fewer 
components,  which  makes  the  alignment  easier.  Despite  the  absence  of  pinholes 
and  saturable  absorbers,  the  ratio  of  ASE  energy  to  pulse  energy  after  the  last 
stage  is  only  a  few  percent.  An  additional  advantage  of  our  system  over  other 
short  pulse  amplifiers,  is  the  fact  that  the  unamplified  dye  laser  pulses  are 
transmitted  through  the  system,  thus  providing  a  c.w.  pilot  beam.  This  is  an 
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important  feature  for  applications  requiring  a  critical  alignment.  Our  ampli¬ 
fier  output  has  a  uniform  spatial  beam  profile,  and  distortion  free  wavefronts. 
The  temporal  properties  of  the  amplified  pulses  have  also  been  investigated, 
usi ng  a  computeri zed  data  acquisition  system,  specially  developed  for  this 
purpose.  The  pulses  are  Gaussian  (FWHM:  6.2  ps),  and  present  a  nonlinear  chirp 
(1.4  E^).  The  bandwidth  is  79  GHz.  It  was  shown  that  the  origin  of  the  chirp 
is  in  the  dye  laser  oscillator,  rather  than  in  the  amplifier  chain. 
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FIGURE  CAPTIONS 


Figure  1.  Schematical  representation  of  the  three  stage  amplifier 
(DBS  =  dichroic  beam  splitter,  POL  =  polarizing  beam  splitter, 
P  »  direct  vision  prism). 

Figure  2.  Spatial  profile  of  the  amplified  beam. 

Figure  3.  Data  acquisition  system:  hardware. 


Figure  4.  Example  of  an  average  fringe,  representing  a  set  of  fringes. 


Figure  5.  Experimental  results  and  best  theoretical  fits  for  (a)  the  intensity 
autocorrel ation  curve  and  (b)  the  fring  resolved  autocorrel ation  curve 
(dashed  lines  are  experimental  curves,  solid  lines  represent  the 
theoretical  fits;  for  the  dotted  line  in  (b),  see  text). 

Figure  6.  Interferometric  autocorrel  ation  curve  for  the  unamplified  pulses  (the 
solid  line  represents  the  best  fit;  the  dashed  line  corresponds  to  the 
best  fit  for  the  amplified  pulses). 
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MULTIPHOTON  IONIZATION  UNDER  SEQUENTIAL  EXCITATION 
BY  COHERENT  PULSES 
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Abstract 

The  use  of  pulse  sequences  offers  a  solution  to  the  conflicting  require¬ 
ments  of  spectral  selectivity  and  high  peak  intensities  needed  to  study  multi¬ 
photon  processes.  In  addition,  it  is  shown  that  random— rather  than  regular- 
pulse  sequences  can  provide  unique  information  on  the  influence  of  coherence  on 
a  resonant  multiphoton  process,  because  the  coherence  can  be  adjusted  independ¬ 
ently  of  the  resonance  condition. 
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I .  INTRODUCTION 


Attention  has  been  drawn  in  recent  years  to  the  possible  use  of  pulse 
sequences  to  enhance  quantum  beat  phenomena  in  molecules,*  c*r  to  selectively 
enhance  weak  spectral  lines. ^  Phase  shifts  in  the  light  field  for  a  single  pair 
of  pulses  can  be  used  to  distinguish  between  single  and  multi  photon  processes.'* 
Statistically  distributed  phase  switching  on  a  CW  field  can  give  it  the  spectral 
and  driving  characteristics  of  a  line  broadened  source.4  Here  we  discuss  an 
application  of  pulse  sequences  in  the  area  of  laser  driven  reactions.  We  choose 
as  an  example  multiphoton  ionization  of  atoms.  However,  the  approach  used  here 
applies  also  to  multi  photon  dissociation  of  mol  ecul es--which  has  been  the  sub¬ 
ject  of  extensive  theoretical  investigation^ — or  to  any  other  process  dependent 
on  the  absorption  of  several  photons. 

Two  types  of  sequential  excitation  will  be  considered:  deterministic  or 
stati stical .  The  motivation  for  using  a  sequence  of  identical  pul ses  with  a 
well  defined  phase  relationship,  is  to  combine  the  spectral  selectivity  offered 
by  a  regular  pulse  train,  with  the  high  peak  intensities  of  picosecond  pulses, 
which  are  needed  to  excite  higher  order  processes. 

A  dimension  of  nonlinear  mechanisms  that  is  subject  to  experimental  test  is 
the  coherence  or  lack  of  coherence  in  the  process.  But,  as  we  discuss  below, 
the  single  pulse  experiments  that  have  been  the  mainstay  of  this  area  provide  a 
fairly  limited  test  of  coherence.  Multimode  versus  single  mode  operation  has 
been  used®’^  to  vary  the  "coherence"  of  the  laser  source.  This  procedure  has 
the  disadvantage  of  varying  the  bandwidth--hence  the  resonance  condition-- 
simul  taneously  with  the  pulse  coherence.  We  will  show  how  a  sequence  of  pulses 
of  controlled  phase  can  be  used  to  vary  the  source  coherence  without  affecting 
its  central  frequency  or  its  spectral  width. 


By  choosing  the  delays  and  phase  between  pulses  to  fit  a  chosen  statistical 
distribution,  we  can  study  how  a  multiphoton  process  is  affected  by  that  statis¬ 
tics.  Such  an  approach,  of  averaging  many  measurements  with  a  predetermi ned 
setting  of  delays  and  phases  (to  fit  a  given  photon  statistics),  will  be  shown 
to  have  several  advantages  over  other  methods  of  studying  the  influence  of 
source  coherence  on  a  multiphoton  resonant  process.  As  an  illustrative  example 
for  the  latter,  we  have  chosen  multiphoton  ionization  in  a  metal  vapor,  in 
presence  of  a  two  photon  resonance. 

II.  THEORETICAL  MODEL 

We  consider  an  atomic  system  in  conditions  of  two  photon  resonance  with 
electromagnetic  radiation  at  u>,  which  can  be  ionized  with  a  single  photon 
transition  from  the  upper  state  (Fig.  1).  The  evolution  equations  for  the 
components  of  the  density  matrix  for  an  exciting  electric  field  of  slowly 
varying  amplitude  can  be  written  (1): 
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T2  is  the  phase  relaxation  time,  ^2  =  2y£2  the  photoionization  rate  from  the 
upper  level  2,  Aw  =  2w  -  -  5u)21  1s  ^recluency  detuning  from  resonance, 

including  the  Stark  shift  <$w^2  With  the  substitution  ia^2  =  Q,  (o22  - 

all)/2  =  W,  and  E2  =  2£2,  and  neglecting  relaxation  and  ionization,  the  equa¬ 
tions  can  be  rewritten: 


Q  =  -iAu)Q  -  W  E2 
W  =  Re(E*  Q) 
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with  the  initial  condition  W(t*0)  =  -1. 


The  equation  (4)  can  be  formally  integrated,  yielding 
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After  a  pulse  has  passed,  Q(  t  >  >  t  )  =  E2p(^),  where  E  2p(  ^ )  is  f-he  Fourier 
transform  of  E^(t).  From  the  conservation  relation  | Q | 2  +  =  1; 


l0l2=  1  -  W2  =  -  2Wqaw  =  2AW,  7 

where  AW  =  W(t  =  “)  -  W(t  =  0).  Combining  Eqs.  (6)  and  (7),  gives  the  energy 
absorbed  by  the  two-level  system: 


2AW  =  -  |E2F(Au>)|2  8 

The  number  of  ions  produced  is  simply  the  number  of  "lost  atoms  **  1  -  °22 
(t  *“>)  -  crj j ( t  =  °°)  in  the  system  of  equations  (l)-(3).  The  number  of  ions  can 
also  be  estimated  by  integrating  the  ionization  rate: 

"ions  “/T lEl%2dt  9 

For  small  population  transfers,  ct22  =  AW,  and 

Nions  =y*v|E|2  A  W  ( t)  dt  10 

Replacing  aW  by  its  average  value  W(°°)/2: 


Nions 
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where  e  =/  |  E|  ^dt  is  the  "pulse  area"  or  tipping  angle  of  the  pseudopol ariza- 

•Coo 

tion  vector  in  the  vector  model  of  the  two-photon  resonant  interaction  (2,3). 

The  expressions  above  hold  in  the  approximation  of  weak  excitation 
(e  small),  no  Stark  shift  (or  any  Stark  shift  but  square  pulse  excitation).  We 
see  that  both  the  absorption  and  the  photoionization  reflect  the  excitation 
spectrum.  The  photoionization  spectrum  is  simply  a  fraction  (proportional  to 
the  "pulse  area")  of  the  Fourier  spectrum  of  the  Rabi  frequency  (which  is,  in 
this  case,  proportional  to  the  square  of  the  exciting  electric  field). 
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III.  COHERENCE  AND  SPECTRAL  CHARACTERISTICS 

We  will  consider  now  more  specifically  the  spectral  characteri sties  (i.e., 
number  of  ions  versus  frequency)  of  an  excitation  mode  of  a  sequence  of  pulses. 
In  general,  even  if  each  individual  pulse  is  pure  in  frequency  (real  amplitude 
function)  the  sequence  of  various  pulses  will  be  represented  by  a  complex 
function  E  (because  of  the  relative  phase  between  pulses).  A  general  property 
of  the  system  of  equations  (4,5)  is  that  W--which  determines  in  first  approxima¬ 
tion  the  photoionization  remains  unchanged  under  a  transformation  E  -*>  E*,  Atu^-Aw 
(indeed,  under  this  transformati on ,  Eq.  4  is  transformed  in  its  complex  con¬ 
jugate  and  Eq.  5  remains  unchanged).  This  symmetry  implies  that  the  photo- 
ionization  spectrum  of  a  single  unchirped  pulse  to  (E  real)  is  symmetric.  The 
symmetry  is  broken  however  in  the  case  of  a  sequence  of  pulses  with  a  phase 
difference  A<>  different  from  0  or  ir,  for  which  we  no  longer  have  E  =  E*. 

It  was  shown  in  the  preceding  section  that  at  least  part  of  the  response  of 
the  medium  may  be  inherent  in  a  Fourier  decomposition  of  the  pulse  sequence. 
Since  we  will  be  considering  sequences  of  pulses  that  are  identical  except  for 
their  relative  phases  and  delays,  the  overall  time  dependence  of  the  signal  can 
be  written  as: 

N  2i*„ 

F(t)  =  £  e  n  (t  +  Tn)  12 
n=l 

In  Eq.  (12)  F(t)  denotes  the  generalized  (complex)  Rabi  frequency  driving  the 
system  of  equation  4,5,  and  f(t)  is  the  Rabi  frequency  (E2)  for  each  individual 
pulses.  The  nth  pulse  is  centered  at  tn,  and  4>n  is  its  absolute  phase.  The 
spectrum  of  F(t)  is  the  function: 

F(n)  =  f(n)  g(n)  13 
which  is  equal  to  the  spectrum  of  a  single  pulse  f(ft)  "modulated"  by  the 
function  g(n)  given  by 


A  particular  case  of  interest  is  the  coherent  sequence  of  N  identical  pulses 
with  relative  time  delay  At  and  (relative)  dephasing  ($n  =  nA$ ;  Tn  =  m), 
which  leads  to  the  particularly  simple  form  for  gU’): 


sin^NAij)  -  4^1  ^ 

sin^$  -  ) 
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The  spectral  intensity  |g(fi)p,  which  was  shown  to  be  in  a  limiting  case  propor¬ 
tional  to  the  ioni  zation  spectrum,  is  ill  ustrated  in  Fig.  2  for  a  N=10  pul  se 
sequence.  The  multiphoton  ionization  from  exactly  this  type  of  pulse  sequence 
will  be  considered  in  the  following  section. 


IV.  IONIZATION  BY  A  REGULAR  PULSE  SEQUENCE 


We  consider  now  the  multiphoton  ionization  produced  by  a  sequence  of  iden¬ 
tical  pulses,  regularly  spaced  in  phase  and  time.  The  system  of  equations  (1-3) 
is  solved  numerically  for  the  four  initial  conditions  corresponding  to  each  of 
the  elements  of  the  density  matrix  set  successively  equal  to  1.  The  matrix 
formed  by  the  four  sets  of  solutions  forms  an  operator  which,  multiplied  by  any 
set  of  initial  conditions,  will  give  the  density  matrix  as  transformed  by  one 
pulse.  Between  pulses,  the  components  of  the  density  matrix  relax  according  to 
the  resonance  condition  and  the  relaxation  times.  The  phase  of  each  successive 
pulse  is  taken  into  account  by  a  simple  transformation  (rotation)  of  the  opera¬ 
tor  matrix.  The  value  taken  by  1  -  -  °22  a^ter  passages  of  the  pulses  is 

the  fraction  ionized. 

Fig.  3  shows  the  Ionization  spectrum  for  a  series  of  10  square  pulses  of 
area  1,  separated  by  100  psec.  The  pulse  duration  is  10  psec.  The  two  spectra 

correspond  to  a  phase  difference  between  pulses  of  -0.5  r  ( - )  or  -.lit 

( - ).  Even  though  the  condition  e<<l  does  clearly  not  apply  for  the  pulse 
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trains,  and  not  even  for  the  individual  pulses,  the  spectrum  qualitatively 
resembles  that  of  the  excitation  pulses.  In  fact,  the  photoionization  curves 
show  all  the  characteristic  features  of  a  Fourier  spectrum:  the  overall  width 
is  related  to  1/t,  an  increase  of  the  phase  difference  between  pulses  shifts 
correspondingly  the  photoionization  spectrum,  and  the  width  of  the  individual 
peaks  in  the  photoionization  spectrum  decreases  with  increasing  delay  between 
pul ses. 

The  Stark  shift  introduces  an  additional  complexity,  by  essentially  cre¬ 
ating  a  different  resonance  condition  during  and  between  pulses.  This  effect  is 
demonstrated  by  the  simple  example  of  a  sequence  of  two  square  pulses  shown  in 
Fig.  4.  Increments  of  delay  between  the  pulses  by  multiples  of  (2ir)/(Stark 
shift),  produce  the  same  number  of  ions  at  resonance  and  at  a  detuning  from 
resonance  equal  to  the  Stark  shift.  Departure  from  the  simple  approximations 
(6-10)  will  be  particularly  significant  for  large  area  pulses.  “Photoionization 
spectroscopy"  with  a  sequence  of  1 arge  area  pulses  will  enable  a  direct  measure¬ 
ment  of  the  photoionization  rate  from  the  upper  state  through  the  width  of  the 
individual  peaks.  Indeed,  as  shown  in  Fig.  5,  the  width  of  each  peak  in  the 
photoionization  spectrum  is  seen  to  increase  with  the  transition  rate  to  the 
continuum.  It  is  interesting  to  note  (Fig.  5)  that  the  peak  ionization  is 
smal  1  er  for  the  higher  ionization  rate.  The  reduction  in  peak  ionization  rate 
is  due  to  the  increased  dephasing  associated  with  the  transition  rate  to  the 
continuum,  resulting  in  a  decrease  in  coherent  excitation  of  the  upper  level. 

V.  IONIZATION  BY  A  RANDOM  PULSE  SEQUENCE 

An  incoherent  field  can  be  described  by  a  density  matrix  corresponding  to 
an  incoherent  superposition  of  coherent  states.®  We  choose  to  use  for  this 
superposition  a  series  of  identical  pulses  with  a  random  distribution  of  delays 
and  phases. 
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Under  certain  conditions,  a  large  number  of  experiments  performed  with  a 
few  pulses,  for  which  the  relative  delays  are  distributed  accordi ng--for  in¬ 
stance— to  a  Poisson  Statistics,  will  have  the  same  average  result  as  the  same 
experiment  performed  once  with  a  very  large  number  of  pulses,  distributed  ac¬ 
cording  to  the  same  statistics.  The  average  time  between  pulses  has  to  be 
shorter  than  the  coherence  time  of  the  particular  resonance  under  investigation. 
We  demonstrate  next  that,  if  the  pulses  are  randomly  distributed,  the  time 
evolution  of  the  system  driven  by  a  sequence  of  identical  pulses  is  proportional 
to  the  number  of  pulses,  and  to  the  total  fluence  delivered,  as  expected  for 
incoherent  excitation.  If,  instead,  the  pul ses  are  appl ied  in  phase  and,  for 
instance,  in  condition  of  two  photon  resonance,  the  absorption  by  the  system  is 
proportional  to  the  square  of  the  number  of  pulses.  For  demonstration  purposes, 
let  us  consider  near  two-photon  resonance  (E2>>Aw)  and  weak  pulses  such  that  the 
approximation  W  =  WQ  can  be  made  in  Eq.  4.  Each  of  the  pulses  of  "area 

.-t  f 

A0  -  I  E2dt  produces  an  increase  in  polarization 
tj 

Q  =  WQ  A0  exp(  i<|> )  16 
If  the  pul  ses  are  al  1  in  phase  (pul  ses  equal  ly  spaced  and  in  phase),  as  is  a 
constant  angle,  and  after  N  pulses,  Q  «  Nao  and  the  energy  absorbed  W  is 
proportional  to  the  square  of  the  number  of  pulses: 


If  the  distribution  of  pulse  delays  is  random  (Awt  d>>l),  then  the  phase  in  Eq.  1 
is  a  random  variable,  and  the  vector  Q  makes  a  random  walk  of  N  steps  aq  in  the 
complex  plane,  starting  from  the  origin.  The  resultant  polarization  is 

Q  *  NAS,  and  the  energy  absorbed  is  proportional  to  the  number  of  pulses: 

A0  2 

W  =  N  18 
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It  should  be  noted  that  by  using  a  sequence  of  pulses,  we  can  compare  coherent 

(as  in  Eq.  17)  versus  incoherent  (as  in  Eq.  18)  excitation,  under  the  same 

conditions  of  total  f  1  uence  and  peak  i ntensi ty.  As  expected,  it  is  only  in  the 

case  of  incoherent  excitation  that  the  evolution  of  the  system  can  be  simply 

proportional  to  the  total  fl  uence  delivered.  The  intensity  dependence  of  the 

multiphoton  process  is  not  invoked  in  comparing  the  incoherent  versus  coherent 

case.  The  above  considerations  sti  1  1  hoi  d  in  the  more  general  case  of  a  n- 

photon  transition.  The  only  difference  is  a  generalized  Rabi  frequency  in  En. 

r  tf 

Equations  16  through  18  remain  unchanged,  with  ao  = /  Endt.  In  the  approxi- 
mation  En>>  Au>,  and  no  intermediate  level  resonance,  coherent  excitation  is 
characterized  by  an  absorbed  energy  proportional  to  the  square  of  the  number  of 
pulses,  while  the  absorption  is  directly  proportional  to  the  number  of  pulses  in 
the  incoherent  case. 

In  order  to  study  the  influence  of  source  coherence  on  the  photoionization 
spectrum  of  lithium,  we  ca  1  cul  a  ted  the  number  of  ions  as  a  function  of  pulse 
wavelength  for  a  coherent  and  an  incoherent  sequence  of  10  pulses  (Fig.  6).  The 
"incoherent"  curve  (solid  line.  Fig.  6)  corresponds  to  an  average  of  500  pulse 
sequences,  each  with  a  different  set  of  relative  delays  taken  from  a  Poisson 
distribution.  The  relative  phases  are  uniform  random  numbers  between  0  and  2n. 
The  total  area  in  each  sequence  of  pulses  is  1.  As  can  be  seen  from  Fig.  6,  the 
photoionization  from  a  coherent  sequence  (10  pulses  equally  spaced,  relative 
phase  ir/2,  dashed  line  in  Fig.  6)  can  be  much  larger  than  that  of  the  incoherent 
excitation.  This  resul  t  may  appear  to  be  in  contradiction  with  experimental 
resul  ts  of  Lecompte,  et  al.  [6],  on  the  influence  of  coherence  on  mul  ti  photon 
excitation.  In  the  case  of  nonresonant  N  photon  ionization  of  atoms  reported  in 
ref  6,  the  ionization  probability  is  proportional  to  the  Nth  power  of  the  laser 
intensity,  and  therefore  very  sensitive  to  the  transient  intensity  peaks  that 


occur  in  a  multimode  laser  pulse.  The  dependence  of  multiphoton  ionization  on 
laser  coherence  was  studied  by  varying  the  number  of  oscillating  laser  cavity 
modes  between  one  and  one  hundred.  It  has  been  estimated  [9]  that,  in  *hat 
case,  a  20  mode  laser  is  a  good  approximation  of  an  incoherent  source.  The 
situation  is  more  complex  in  the  case  of  resonant  multiphoton  ionization  as  we 
are  considering  here.  Indeed,  the  number  of  osci  1  1 ating  modes  i s  varied  by 
changing  the  bandwidth  of  the  cavity.  Therefore,  the  parameter  "coherence" 
cannot  be  controlled  independently  of  the  resonance  condition,  since  the  spec¬ 
tral  width  of  the  laser  source  increases  with  decreasing  coherence.  In  the 
approach  used  here  however,  the  laser  coherence  is  modified  at  constant 
bandwi dth,  since  the  spectrum  of  the  pulse  sequence  is  contained  within  the 
spectrum  of  a  single  pulse. 

It  is  interesting  to  note  that,  in  a  recent  experimental  study  of  4  photon 
(3  photon  resonant)  ionization  of  caesium  atoms,  Lompre,  et  al.  [7)  found  also 
an  increase  in  the  width  of  the  spectral  dependence  of  the  ionization,  for 
incoherent  versus  coherent  excitation.  The  hyperfine  structure  of  the  resonance 
is  still  observed  in  coherent  excitation,  while  it  is  not  in  the  incoherent  case 
(ref.  7,  Fig.  8).  However,  in  that  particular  experiment  [7],  the  coherence  was 
varied  by  changing  the  number  of  modes  of  the  laser,  and  comparing  data  for  the 
same  average  intensity.  Therefore,  temporal  fluctuations  make  the  peak  i nten- 
si_ty  highest  in  the  incoherent  case,  resulting  in  a  higher  ionization  in  the 
incoherent  case  than  in  the  coherent  case.  In  an  experiment  in  which  the 
coherence  is  varied  by  varying  the  temporal  distribution  of  identical  pulses,  as 
suggested  here,  the  peak  intensity  as  well  as  the  spectral  bandwidth  remain 
constant,  while  the  coherence  is  modified  within  the  pulse  spectrum. 


Photon  Statistics 

Let  us  consider  a  limiting  case  where  each  pulse  of  the  sequence  contains 
only  one  photon.  The  "pulse  sequence  statistics"  is  then  nothing  else  than  the 
photon  statistics  of  the  radiation  field.  The  number  of  pulses  N  is  now  the 
number  of  photons  in  the  excitation,  each  of  these  photon  -h™  driving  coherently 
a  n-photon  resonant  system  (i.e.,  a  two  level  system  of  energy  n -hoi).  For 
instance,  in  the  simple  argument  at  the  beginning  of  this  section,  the  elemen¬ 
tary  increment  of  the  n-photon  off  diagonal  element  becomes 


A9  - J  Endt  =  [-ho)]n/2  19 

n  being  a  proportionality  coefficient. 

It  should  be  noted  that  the  photon  energy  does  not  appear  explicitly  in  any 
calculations  carried  out  in  the  slowly  varying  envelope  approximation  (SVEA). 
The  statistical  sequence  of  identical  macroscopic  coherent  pulses  (each  of  them 
being  a  coherent  superposition  of  M  photons)  interacting  with  the  n-photon 
resonant  atoms,  is  identical  (within  the  SVEA)  to  a  photon  statistics  in  a 
ficti  tious  worl  d  where  PI  anck's  constant  is  M-h-.  The  transition  frequency  at 
u)0  =  no>  is  in  condition  of  n-photon  resonance  with  each  pulse  (Mb)<»>0  =  n(Mt»u>). 
Since  Planck's  constant  appears  only  implicitly  (in  the  scaling  of  the  gen¬ 
eralized  Rabi  frequency  E)  in  the  theoretical  model,  all  the  equations  of  the 
theoretical  model  remain  unchanged  by  the  "transformation."  The  elementary 
increment  of  the  off  diagonal  el  ement  for  the  standard  pul  se  of  the  train  can 
also  be  expressed  in  terms  of  the  "pseudo  photon"  energy  (Mh)oj,  in  a  similar  way 
as  in  Eq.  19: 


A0  ~J" Endt  =  n[(Mh-)u>]n/2 


20 


This  similitude  offers  an  intuitive  justification  to  the  use  of  a  pulse  sequence 


to  study  the  statistical  response  of  a  multiphoton  resonant  system.  We  are 


esen t i ally  creating  a  pseudo  photon  statistics.  This  concept  is  legitimized 
because  PI  anck's  constant  does  not  appear  expl  icitly  in  this  formal  i  sm  ( i.e., 
the  photon  energy  i s  sti 1  1  small  compared  wi  th  the  tota  1  exci tation  energy). 
Since  the  model  of  the  pseudostatistics  concept  is  equivalent  to  the  theoretical 
model  used  in  section  III,  it  does  not  bring  any  new  light  to  the  problem. 
However,  it  offers  a  convenient  i  nterpretation  to  the  approach  that  is  being 
taken  here,  to  study  coherence  effects  through  the  use  of  pulse  sequences. 

The  work  was  supported  by  the  Air  Force  Office  of  Scientific  Research 
(AFOSR)  under  Grant  No.  820332,  and  by  the  Robert  A.  Welch  Foundation. 
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FIGURE  CAPTIONS 


Figure  1.  Level  diagram. 

Figure  2.  Plot  of  the  spectral  intensity  for  a  sequence  of  10  identical  square 
pul  ses,  as  defined  in  Eq.  15. 

Figure  3.  Ionization  versus  angular  frequency  (ns-*)  for  a  sequence  of  10  Gaus¬ 
sian  pulses,  with  an  interpulse  separation  of  0.1  ns.  The  area 
(energy)  of  each  indi  v idua  1  pu  1  se  is  1 ,  and  thei r  duration  ( FWHM  of 
the  intensity)  is  0.01  ns.  Two  different  cases  of  relative  phase 

difference  are  shown:  =  -  ir/2  ( - )  and  =  -0.1  ( - ).  The 

parameters  are  taken  from  the  2s-4s  transition  of  lithium  (Stark 
shifts  -  1.5  E2:  ionization  rate  from  the  4s  level:  0.07  E2). 

Figure  4.  Ionization  versus  angular  frequency  (ns-*)  for  a  sequence  of  two 
square  pulses.  The  individual  pulses  have  a  duration  of  0.05  ns  and 
an  amplitude  of  E  =  10  ns"*/2.  The  medium  parameter  are  the  same  as 
for  Fig.  2  (hence  the  Stark  shift  is  -150  ns"1).  The  solid  line  shows 
the  ionization  spectrum  for  a  sequence  of  two  pulses  90°  out  of  phase 
following  each  other  without  delay.  For  the  dashed  line,  the  two 
successive  pulses,  90°  out  of  phase,  are  separated  by  a  delay  of 
[2  /(Stark  shift)]. 


Figure  5.  Single  peak  of  the  ionization  spectrum  (%  ions  versus  angular  fre¬ 
quency  in  ns"*),  for  two  different  ionization  rates.  The  exciting 
field  consists  in  a  sequence  of  10  Gaussian  pul  ses,  of  0.01  ns  dura¬ 
tion  (FWHM),  and  area  e  *^E2dt  *  3.  The  successive  pulses  are  0.4 
out  of  phase.  The  medium  parameters  are  the  same  as  in  the  preceding 
figures,  except  for  the  Ionization  rates  YE2  being  considered  here. 
The  scale  to  the  right  corresponds  to  the  ionization  rate  of  0.07  E2 


( - )  while  the  scale  to  the  left  pertains  to  the  higher  ionization 

rate  of  0.44  E^  ( - ). 

Figure  6.  Ionization  versus  angular  frequency  (ns-1)  for  a  sequence  of  10 
pulses.  Each  pulse  is  Gaussian  shaped  of  area  0.1.  The  solid  line 
corresponds  to  a  random  pulse  sequence;  the  dashed  line  to  a  sequence 
of  equally  spaced  (0.09  ns)  pulses,  each  90°  out  of  phase  with  the 
former.  Each  pulse  of  the  sequence  is  Gaussian,  with  a  duration  of 

0.01  ns  (FWHM).  The  statistical  result  ( -  )  is  obtained  by  taking 

the  average  of  500  "excitations,"  each  with  a  sequence  of  pulses 
distributed  (in  relative  delay)  according  to  a  Poisson  distribution. 
The  relative  phases  between  pulses  are  taken  to  be  random  numbers 
between  0  and  2  .  The  medium  parameters  are  the  same  as  for  Fig.  2. 
The  average  delay  between  pulses  is  0.09. 
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Figure  3. 


IONIZATION  FOR  A  10  PULSE  SEQUENCE 


IONIZATION  FOR  A  2  PULSE  SEQUENCE 
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t  DELAY 


FREQUENCY  FROM  RESONANCE 
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Figure  4. 


Sequence  of  10  Pulses  90  ’  out  of  Phase 
Random  Sequence  of  10  Pulses 


INTERFEROMETRIC  AUTOCORRELATIONS  APPLIED  TO  THE  STUDY  OF 
MULTI PHOTON  COHERENCE 


H.  Vanherzeele  and  J.-C.  Diels 
Center  for  Applied  Quantum  Electronics 
Department  of  Physics 
North  Texas  State  University 
P.O.  Box  5368;  NT  Station 
Denton,  Texas  76203 
817-565-3308 


A  new  method  to  study  multiphoton  coherent  excitation  by  intense  pico¬ 
second  pulses  is  presented.  The  experimental  setup  includes  a  dye  laser 
oscillator-amplifier,  an  interferometric  autoco^rel  at.or  and  a  computerized 
data  acquisition  system. 


INTERFEROMETRIC  AUTOCORRELATIONS  APPLIED  TO  THE  STUDY  OF 
MULT I PHOTON  COHERENCE 


H.  Vanherzeele  and  J.-C.  Diels 
Center  for  Applied  Quantum  Electronics 
Department  of  Physics 
North  Texas  State  University 
P.0.  Box  536R;  NT  Station 
Denton,  Texas  76203 
R1 7-565-3308 

We  developed  a  high  gain  picosecond  pulse  dye  laser  osci  1 1 ator-ampl  ifier 
system,  coupled  to  an  interferometer  and  a  computerized  data  acquisition 
system,  to  perform  nonlinear  experiments  with  accurately  phased  pulse 
sequences  in  vapors.*  Data  acquisition  and  processing  is  performed  on  line  at 
the  10  Hz  repetition  rate  of  the  laser  source.  The  interferometric  delay  line 
splits  the  main  pulse  into  two  pairs  of  pulses  with  continuously  increasing 
delay.  The  relative  phase  between  the  pulses  (slowly  increasing  at  a  rate  of 
circa  2t>/4s),  as  well  as  the  absolute  delay  between  the  pulses  is  determined 
from  a  HeNe  laser  i nterferogram.  Our  setup,  which  permits  to  analyze  the 
waveform  and  phase  modulation  of  the  pulse,  is  designed  to  investigate  multi- 
photon  coherent  effects  by  cross  correlating  the  mul ti photon  stimulated  rera¬ 
diation  with  the  pulse  itself. 

The  ratio  of  the  second  harmonic  after  the  interferometer  to  that  before 
the  interferometer,  provides  both  the  intensity  and  interferometric  autocor¬ 
relation.  These  measurements  enabled  us  to  determine  the  pulse  shu^c  tr  be 
E  ■  exp  r -( t/5.27)2  ♦  i 01  with  0  =  1.43  E^,  and  t.  in  ps. 

The  other  pair  of  pulses  is  sent  in  conditions  of  multiphoton  resonance 
into  a  vapor.  The  second  harmonic  of  the  transmitted  signal  is  recorded  as  a 
function  of  the  delay  and  the  relative  phase  between  the  pulses.  The  perio¬ 
dicity  in  phase  of  the  transmitted  signal  is  a  direct  measurement  of  the  order 
of  the  process  (2nfor  single  photon,  n  for  two  photon,  V2  for  four  photon 


Vanherzeele,  et.  al . 

Interferometric  Autocorrel ations  Applied  to... 

resonance).  The  shape  and  decay  of  the  signal  provides  simultaneously  infor¬ 
mation  about  coherence  times  as  well  as  spectroscopic  information. ^  It  should 
be  noted  that,  because  of  our  accurate  knowledge  of  the  pulse  shape,  our 
temporal  resolution  is  at  least  one  order  of  magnitude  smaller  than  the  pulse 
duration  (6.2  ps  FWHM). 

The  software  of  our  data  acquisition  system  averages  out  data  over 
several  periods  for  a  given  phase,  in  order  to  generate  an  accurate  display  of 
signal  (or  transmission)  versus  phase.  Simultaneously,  phase  averaged  values 
of  the  data  are  also  computed.  As  an  example  of  the  latter,  the  second 
harmonic  of  the  signal  transmitted  through  Li  vapor,  as  a  function  of  delay 
between  the  pulses,  is  displayed  in  Fig.  1.  The  laser  pulses  are  tuned  to 
resonance  with  the  2s-4s  two  photon  transition  of  Li  at  571.2  nm.  In  the 
absence  of  Li,  the  averaged  data  correspond  to  the  conventional  3:1  intensity 
autocorrel ation  curve.  The  enhanced  averaged  attenuation  in  Li  near  zero 
delay  is  in  perfect  agreement  wi th  a  quantum  mechanical  calculation  of  the 
coherent  interaction  of  the  radiation  with  the  two  level  system.  The  phase 
dependent  data  on  the  other  hand,  are  much  more  complex  but  also  richer  in 
information.  Examples  of  coherent  interaction  in  dye  vapors  (single  photon), 
Li  vapor  (two  photons),  and  Hg  vapor  (four  photons)  will  be  presented  to 
demonstrate  this. 
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